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We re- visit the production of cosmic rays by cusps on cosmic strings. If a scalar field ("Higgs") has 
a linear interaction with the string world-sheet, such as would occur if there is a bosonic condensate 
on the string, cusps on string loops emit narrow beams of very high energy Higgses which then 
decay to give a flux of ultra high energy cosmic rays. The ultra-high energy flux and the gamma to 
proton ratio agree with observations if the string scale is ~ 10 13 GeV. The diffuse gamma ray and 
proton fluxes are well below current bounds. Strings that are lighter and have linear interactions 
with scalars produce an excess of direct and diffuse cosmic rays and are ruled out by observations, 
while heavier strings (~ 10 15 GeV) are constrained by their gravitational signatures. This leaves a 
narrow window of parameter space for the existence of cosmic strings with bosonic condensates. 



The detection of cosmic topological defects can pro- 
vide a direct window to fundamental physics and the very 
early universe (for a review, see [J). Hence it is no sur- 
prise that there has been a concerted effort to examine 
observational consequences of cosmic topological defects, 
including gravitational wave emission, gravitational lens- 
ing, fluctuations in the cosmic microwave background, 
and ultra-high energy cosmic rays. If cosmic strings are 
superconducting, they may also lead to electromagnetic 
phenomena such as gamma ray and radio bursts. 

The gravitational effects of cosmic strings are stronger 
for heavier strings and current observations of the cos- 
mic microwave background rule out strings above the 
~ few x 10 15 GeV energy scale. As lighter strings are 
considered, gravitational effects become less significant, 
and other particle physics signatures become relatively 
important. If the strings are superconducting, the cur- 
rents on the string can lead to electromagnetic radiation 
that could be observable. However, if the strings are not 
superconducting, such signatures will be absent, and one 
must turn to particle emission from cosmic strings. If 
string dynamics forces loop production at the smallest 
possible scales, as suggested in 0], particles can be co- 
piously emitted, leading to strong constraints. However, 
other studies indicate that string loops are large com- 
pared to microscopic length scales [l[ and particle emis- 
sion is suppressed Q. Even in this case, portions of a 
string loop may get boosted to very high Lorentz factors, 
creating a "cusp" on the string (see Fig. [T|) , and this may 
potentially provide a burst of particles that could be seen 
in cosmic ray detectors. 

Particle emission from cusps on cosmic strings has been 
studied by several authors. Srednicki and Thcisen Q 
considered a quadratic interaction of a scalar field with 
an idealized (zero thickness) string and came to the con- 
clusion that particle emission is insignificant for astro- 
physical size string loops. Our analysis for the radiation 
is similar to that of Refs. [J, [f|, though the particular 
linear interaction of the scalar field (call it H) with the 
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FIG. 1: String segment at 3 different times, with a cusp at 
O, with velocity along z. The size of the curved section of 
string is L. The observer is at an angle 8o from the z-axis 
at a distance d. There is strong emission from the cusp in a 
beam along the z axis. 



string world-sheet, as would occur if H condenses on the 
string (27|, has not been considered before. A linear in- 
teraction causes an enhancement of particle production 
by a factor of M/m where M is the string scale and m is 
the mass of H. If M is the Grand Unification scale, while 
m is the electroweak scale, this factor can be as large as 
10 13 . 

Particle emission from cusps on thick strings has been 
considered in Ref . [|| 0| ■ Now the cusp consists of over- 
lapping, oppositely oriented strings that can annihilate 
and give off energetic particles. A careful study of this 
process, including numerical evolution of the field equa- 
tions, shows that the resulting flux of particles is too 
small to be of interest [!,[§]. In contrast, the linear inter- 
action mechanism wc study is insensitive to the thickness 
of the string, and occurs over a string length that is much 
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larger than the length over which cusp annihilation oc- 
curs. Thus we can ignore cusp annihilation and work in 
the zero thickness limit. 

There are important observational constraints on the 
so-called "top down" models for production of ultra-high 
energy cosmic rays (UHECR), in which a heavy particle 
decays to give ultra-high energy protons and gamma rays. 
The constraint arises because there are bounds on the 
diffuse gamma ray flux in the EGRET window and also 
on the fraction of photons to protons in the observed 
UHECR. In previous studies, using cusp annihilation on 
cosmic strings as the source of Higgs injection, it was 
found that the EGRET constraint on diffuse gamma ray 
fluxes implies that the flux of UHECR is uninterestingly 
low. However, the constraint on diffuse gamma rays is 
sensitive to the spectral features of the injected Higgs 
particles and to the particular interaction of the string. 
The Higgs particle emission that we consider yields a 
diffuse flux of gamma rays and protons that is below 
EGRET bounds if the string scale is > 10 6 GcV and 
an UHECR flux consistent with observations for string 
scale ~ 10 13 GeV. The ratio of photons to protons in 
the UHECR flux is also consistent with current bounds. 
Hence our model can explain the observed UHECR and 
not run into trouble with the EGRET bounds. 

It is to be noted that cosmic ray production by strings 
increases as the string scale decreases. Hence strings on 
scales less than ~ 10 13 GeV, and with these interactions, 
are excluded by the observed flux of UHECR. This is 
especially interesting since heavier strings are constrained 
by their gravitational signatures. Hence there is a narrow 
window between, say, 10 13 — I0 15 GeV for the mass scale 
of cosmic strings having linear interactions with a scalar 
field. 

To summarize, the novelty in the present work is that 
we are considering a new, generic, interaction of cosmic 
strings with scalar fields that leads to high energy parti- 
cle emission. This interaction seems to have been missed 
in the literature. Also, we have been careful to take the 
beamed nature of the emission into account. We have 
focused on deriving approximate analytical estimates so 
as to keep the physical aspects of the problem as appar- 
ent as possible. More detailed predictions will require 
numerical evaluation of the production and propagation. 

We begin in Sec. [J by describing the field theory inter- 
actions under consideration. We then evaluate the rate 
at which a single cusp emits Higgs particles in Sec. [n] In 
Sec. IIIII wc use the results of Sec. [H] to determine the cos- 
mological Higgs injection function. The diffuse gamma 
ray and proton fluxes are calculated in Sees. IIVI and [V] 
respectively. These cosmic rays originate in strings that 
are relatively far away from the Earth. Higgs emission 
from strings that are closer to us and pointed at us can 
give UHECR. The direct flux of UHECR is calculated 
in Sec. I VII where we also discuss the ratio of photons to 
protons in UHECR. Wc conclude in Sec. IVIII In Ap- 
pendix [A] wc summarize some known facts about cosmic 
string loop dynamics. 



I. FIELD THEORY 

The interaction we consider is 

S int = -kM J d 2 ay/^ h (1) 

where k is a coupling constant assumed to be ~ I, M 
is the string energy scale assumed to be at the Grand 
Unified scale, "f a b the string world-sheet metric, and h a 
scalar field. 

A linear interaction can be considered quite generally. 
It can also arise if there is a bosonic condensate on the 
string. To see this explicitly, consider the model 

S = S [$,H,...]+k J d 4 x(&® - M 2 )H^H (2) 

where So is a field theory action that yields cosmic 
string solutions when $ gets a vacuum expectation value 
(VEV). H is a scalar field that we can take to be the 
electroweak Higgs for concreteness. 

At energy scales above the Grand Unified scale, the 
VEVs of $ and H are both zero. At lower energy scales, 
but still above the electroweak scale, $ gets a VEV so 
that |($)| 2 = M 2 but \(H)\ 2 = 0. At this stage, we 
also have string solutions whose tension is fi ~ M 2 and 
width is ~ M _1 . Inside the core of the string, where 
can become small, it may be favorable for H not to 
vanish since the coefficient of the H^H term in Eq. ^ 
becomes negative. As in the case of bosonic supercon- 
ducting strings [lo| . there can be an H condensate in 
the core of the string. Since M is the only mass scale 
in the problem at this stage, the magnitude of H is of 
order M within the core of the string. At a lower en- 
ergy scale, H too gets a VEV. For example, if H is the 
electroweak Higgs, this scale is m ~ 100 GeV. We will as- 
sume m <C M and hence the VEV of H within the string 
is unaffected by the lower scale (electroweak) symmetry 
breaking. 

The interaction term in Eq. ^ can now be written as 
S int = kJ d 2 a J d 2 x±</=7(&$ - M 2 )H^ H 

= k J (fay/^f J d 2 x±(&$-M 2 ) x 

((H) in + h^((H) iD + h) 
« -kM J d 2 aV^h+... (3) 

In the first line we have split the integral into world-sheet 
and transverse integrals and included the Jacobian factor 
(1/— 7) where j a b (a, b = 0, 1) denotes the induced metric 
on the string world-sheet. The subscript "in" on the an- 
gular brackets denotes that the relevant value is the VEV 
within the core of the string and we take (-ff)m ~ M. 
Note that h in the last line denotes the radial component 
of H evaluated on the world-sheet. 

In the next section, we will estimate the radiation of 
Higgses from cusps on cosmic string loops. 
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II. HIGGS EMISSION 



The equation of motion for the Higgs held is 



(D + m 2 )h 



(4) 



where 



3(x) 



kM / dr da 



-7 6 4 (x- X(<j,t)) (5) 



Then the number of Higgs particles with momentum k 
produced due to a source is 



dN k = \j(co k ,k)\ 2 p- 
2uj k 



(0) 



where the Fourier transform of the source j(x) is given 
by 



j(wfe,k) = —kM / drdayf—^e 



i[w fc T-k-X(cr,r)] 



(7) 



where X^(a, t) denotes the string world sheet and ui k = 
\Jk 2 + m 2 . 

The dominant contribution to j comes from the region 
around the cusp where \u k r — k • X(er, r)| < 1. Choosing 
world-sheet coordinates so that the cusp occurs at a = 
and t = 0, this occurs for a range |t|, \a\ < L/(kL) 1 ' 3 
provided k > m\JmL (see Appendix [X]). Also, 



Therefore 



-=l-x 2 ~^L~(fcL)- 2 / 3 



(8) 



(9) 



The angular width of the beam of particles emitted from 
the cusp can also be estimated by evaluating the integral 
in Eq. (f71) in the stationary phase approximation, or as 
in Ref. The result is 



(10) 



(fcL)V3 

Then Eq. ([6]) gives, 

dN k ~ \j\ 2 2 kdk ~ n 2 M 2 L 2/3 -^ 
The estimate applies for 

k G fmV mL, mVml) . 



dk 



k 7/3 



(11) 



(12) 



The upper cut-off on k arises because the wavelength of 
the emitted particles should be larger than the string 
width in the rest frame of the cusp: A > M -1 . Boosting 
to the rest frame of the loop, this yields k < M/y/1 — x 2 
and together with the estimate in Eq. ((5J) gives the upper 
cut-off. The lower cut-off comes from the requirement 




cusp annihilation 



FIG. 2: Higgs emission with momentum k occurs from a sec- 
tion of the string with length ~ Lf(kL) 1 ^ 3 . Cusp annihilation 
occurs over a length ~ y/L/M [T3] and is smaller than the 
region for Higgs emission for k < M V ML. 



that \k-X\ < 1 (see Appendix IA"]). For k < mymL, \a\ ~ 
k/m 2 and the spectrum is a rapidly increasing function 
of k 



dN k 



M 2 k 25 ? 3 

m 16 L 14/3 



dk 



(13) 



Hence dN k goes to zero very fast as k — > 0. In the fol- 
lowing sections, wc will ignore this part of the spectrum 
and only consider k > m\JmL. 

At this stage we can also compare Higgs emission from 
cusps to the process of cusp annihilation (see Fig. 
Higgs emission at momentum k occurs over a length 
Lf(kL) 1 / 3 . With the upper bound, k = /c max = M V ML, 
this length is y/L/M and coincides with the cusp anni- 
hilation length [13] ■ Hence, cusp annihilation docs not 
affect our estimates of Higgs emission for k < fc max . 

A caveat to this statement is that the presence of the 
condensate should not significantly change the string dy- 
namics. Also, we have been considering cusp annihila- 
tion, but the condensate itself has some width which is 
larger than the string width and there could, in principle, 
be "condensate annihilation" even where there is no cusp 
annihilation 1131. 



III. HIGGS INJECTION FUNCTION 

Let d^n(k,t) be the number of Higgs produced by 
string cusps with energy (k,k + dk) at time t per unit 
volume per unit time. d<S>H is called the "Higgs injection 
function" . Then 

/, x „ /* imax ,r dN H dn L dN c . „. 

where the first factor in the integrand is the number of 
Higgses with energy k produced by a cusp on a loop of 
length L. The second factor is the number density of 
loops of length L at time t. The third factor is the num- 
ber of cusps per unit time on a loop of length L. The 
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integration is over all loops of length from L mm to L max - 
The smallest loop can have length ~ M _1 but most of the 
ultra-high energy cosmic ray signal will come from longer 
loops. L max is clearly bounded by the cosmic horizon size 
~ t but, for a fixed value of k, it is also bounded due to 
the constraints in Eq. (fT2|) . 

Let us deal with the last factor in (fl4|) first. Since the 
motion of the loop is periodic or quasi-periodic 



dNc 
dt 



k 

L 



(15) 



where f c is a parameter that gives the average number 
of cusps on a loop per oscillation period. For loops that 
aren't too complicated, we expect f c ~ 1. 

The rate of Higgs production from a single cusp - the 
first factor in the integral in (fT4")) - is given in Eq. (fTTj) . 
Now we need to determine the second factor - the number 
density of loops - in the integrand in Eq. (fT4|) . 

The number distribution of cosmic strin g lo ops is cur- 
rently under discussion @, Q3, GJ, M, El El In the 
authors find that loops will only be of microscopic size, 
in which case loops decay very quickly. The conventional 
scenario, though, is where there is a distribution of loops 
of all sizes at all times. Simulations of (non-radiating) 
string networks in an expanding universe give 



dL 



A 



(16) 



where A ss 10 and Li is the initial length of the loop. 
In these simulations, the loops do not shrink due to ra- 
diation and L = Li at all times. To include the effects 
of radiation from string loops we will express the initial 
length, Li, in terms of the length at time t. This is given 
by the differential equation for the rate of energy loss 



dL „ „ 



/' 



(17) 



where T g and Th arc numerical coefficients characterizing 
the gravitational and Higgs radiation. The energy lost 
to Higgs radiation is found by integrating Eq. (fTTj) after 
multiplication by k. In Eq. (|17p . the energy lost in one 
cusp event is averaged over one oscillation period of the 
loop to get a rate of energy loss. 

Eq. (fT7|) has to be solved with the initial condition 
L{t = 0) = Li, which is equivalent to assuming that 
all the string loops are effectively laid down at t = 0. 
Strictly, there will be loops of size Li that are produced 
at some later cosmic time ti, and these will then evolve. 
However, for the purpose of finding the distribution func- 
tion for loops, one can view these as larger loops that were 
produced at t = which then shrunk to Li at the time 
ti. This is the usual scheme for determining the loop 
distribution function with the inclusion of loop evapo- 
ration (e.g. see [l|). A more sophisticated analysis that 
attempts to take radiation backreaction into account fl8j 
gives a similar loop distribution function though with a 




FIG. 3: The integration in I(Y, x) is over y (dashed vertical 
line) at a definite value of x (say xo). The different loop 
trajectories are labeled by the initial length, yi, at i = 0. 
The dotted line at y = 1 is where the differential equation 
changes behavior (see Eq. (J20J ) . For small values of y(xo), 
yi ~ x up to factors of order 1. For large y(xo), yi ~ y. 



different dependence on L. However, the effects of back- 
reaction on the loop production function, which feeds into 
the loop distribution function, has not yet been included 
[l8| and hence we have stayed with the conventional loop 
distribution function. 
With the rescalings 



y 



Eq. (JTTJ) becomes 
dy 



mL , x 



T h mt 



(18) 



(19) 



The differential equation can be solved but it is simpler 
to approximate it as 



dx 
dx 



= -i , y>i 
i 



, y < i 



(20) 



The Higgs injection function in Eq. (Q3J) can now be 
written as 



t 2 



4/3 7^(») 



where Y = otL r 



with 



T g Gii 



(21) 



(22) 
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and 



Then 



I(Y,x) 



dy 1 

o y 1/3 v? 



(23) 



Note that we have set the minimum length to since 
the integration will be dominated by the upper limit. 
The maximum length (needed to determine Y) cannot 
be more than the Hubble size but it is also restricted by 



the value of k since k > fc n 



Y 



Therefore 



*1 

TTV 3 



(24) 



To determine I(Y, x) we need to write yi in terms of 
y and x in Eq. (|2"3")) . For this we need to solve Eq. (|2T)|) 
for all yi to obtain y(x\ yi). Then this function should be 
inverted to get yi{y{x), x) which can then be inserted in 
Eq. (|2"3"|) (see Fig. [3]). We shall do these steps approxi- 
mately but analytically. 

If Y is very small compared to x, the integration over 
y is for y <C x and then yi at x up to some numerical 
factor of order 1. Therefore 



y2/3 

J(Y, a:) w —5- , r«i 



(25) 



In the opposite limit, Y ^> x, the integral needs to be 
split into two pieces, one for y from to x and the other 
from x to Y. The first piece is approximated as for the 
Y <C x case. The second piece is different since here 
yi ~ y is more appropriate. More explicitly, for Y x, 



I(Y,x) 



=f 


dy 


1 




Jo 


yl/3 






x 2/3 




r 


dy 


X 2 


+ J 


X 


yl/3 


1 








~ x i/3 









(26) 



Note that I(Y, x) does not depend on Y and, since Y is 
k dependent (Eq. (124j) ). neither does I(Y, x) depend on 
k in this limit. 



Therefore 



d$ H {k,t) 



i 2 



v 4 /3 



fc 7/3 



and 



d$ H (k,t) 



M/e 4/3 
i 2 



dfc 1 

F/ 3 x 4/3 



F < a; 



(27) 



(28) 



A more rigorous derivation will lead to a smooth inter- 
polation between these asymptotic forms. For the pur- 
pose of our estimates, it is sufficient to extend the above 
asymptotic forms so that they connect continuously at 
Y = x. 

These equations can be written more neatly in terms 

of 



iy/T gG[imt 



(29) 



d* ff (fc,t) 
d$ H (k,t) 



/iAf c m 4 dk 

t 2 ~kfT ' 

fj,Af c m 4 dk 



k<K (30) 
K > k>K (31) 



t 2 1.8/3 £7/3 ' 

where K = ra\J mt and the bound k < K ensures that 
Y = ak 2 /m 3 i.e. the loops are less than the horizon 
size. However, the k~ 7 ^ 3 fall off is also valid for k > 
K because then I(Y, x) in Eq. (f2"Tj) is independent of k. 
So the restriction K > k in Eq. (f3Tj) may be dropped. 
Eq. (|31[) is our estimate for the Higgs injection function. 

An important feature of the injection function in 
Eq. (|31[) is that it is inversely proportional to the string 
scale. Thus lighter strings inject more Higgses than heav- 
ier strings. This feature was also noted in Ref. Q and 
can be explained by the greater longevity of light string 
loops. 



IV. DIFFUSE PHOTON FLUX 

Once Higgses are injected into the cosmological 
medium, they will decay, lose energy and eventually cas- 
cade into gammas in the EGRET energy range. The 
energy in gammas in the EGRET window is estimated 
fraction of the total injected energy using 



L 

2 



dt / dk k- 



H 



1 



dk (1 



(32) 



where f v is the fraction of energy of the Higgses that goes 
into pions and z is the cosmological redshift. The 1/2 
accounts for how much energy goes into gammas. The 
diffuse gamma ray background measured by EGRET is 
on the order of w C as,obs ~ 10~ 6 eV/cm 3 in the energy 
range 10 MeV-100 GeV 0. 

With the Higgs injection function found in the previous 
section, we have 



dkk 



di> 



H 



dk 



k * dkt iAf c m 4 



t 2 

nAf c m 



i 



(r 5 G^m) 3 / 2 i 7 /2 



(33) 



The integration is dominated by the upper limit (fc»). 
The integration for k <E (&*, oo) will be dominated by its 
lower limit because of the faster (fc~ 4 / 3 ) fall off and will 
give a comparable result. 

Now with 1 + z = a(to)/a(t) where a(t) is the scale fac- 
tor in a matter dominated universe the time integration 
can be done and gives 



3Afirf c fi 1/2 m 
{Y g G^/ 2 m 2{ °> i 3 



(34) 



The time integral is dominated by the present cosmic 
time and Higgs production during the radiation domi- 
nated era has been ignored. As noted at the end of the 



previous section, the energy cascade into gamma rays is 
larger for lighter strings. 

To get a feel for the numbers, we use m = 10 2 77i2 GcV, 
H = (10 13 7?i 3 GeV) 2 (therefore Gfx = lO" 12 ?? 2 ,), to = 10 28 
cm. Also we take A ~ 10, T g ~ 100 and we assume 
U ~ 1. Then 

^cas - lO-^^ma 1/2 (35) 
cm J 

This is quite a bit smaller than the EGRET observation 
unless the string scale is less than ~ 10 6 GeV. Lighter 
strings with the interactions in Eq. $2$ are ruled out, 
though with the caveat that the dynamics of much lighter 
strings (M < 1 TcV) can be dominated by friction until 
the present cosmic epoch and so the network properties 
and loop dynamics can be quite different 
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FIG. 4: A cusp on a loop close to Earth emits particles that 
lose energy as they propagate and also spread out. 



VI. DIRECT FLUX 



V. DIFFUSE PROTON FLUX 

We shall follow Berezinsky ct al [23] and use a power- 
law fragmentation function to obtain the diffuse proton 
flux, I P (E), as (see Eq. (A13) of (H) 



I P (E) 



(2 - p)f N h H 
Airp E 



E 
Eh~ 



R P {E) (36) 



Here p = 1.9 |2J] and Jn is the fraction of energy trans- 
ferred to nuclcons when a Higgs decays which, for order 
of magnitude estimates, we will take to be 1. fiH is the 
rate at which Higgses are being produced per unit vol- 
ume and can be found from our result for d<&n above. 
Eh is the energy at which the Higgses are produced and 
R P (E) is the proton attenuation length at energy E due 
to scattering off the CMB. At energies of 10 19 eV i.e. 



below the GZK cut-off, R p ~ to 
To estimate I P (E) wc take 



10 



28 



E H = fc* 



10 



19 



r?i 3 m 2 



3/2 



GcV 



(37) 



10~ 58 m2 2 ??r3 2 cm ~ 3 s ~ isr 



3 „-!„„-! 



(38) 



and 

This gives 

E 3 I p (E) ps 10 1 

(39) 

where E\g = 10 19 eV. This is to be compared with the 
observed flux 12511 



jjAfc m 4 

L K * 



E 

Eig 



1.1 



-0.65-1.1 
! '/13 



eV 2 



m^ — s — sr 



[Ef 9 I p (E 19 )] ohs PS 10 5 



eV 2 



sr 



(40) 



Hence the proton flux is much less than the observed 
flux unless r/ 13 ~ 10 -6 , in which case both the diffuse 
gamma flux and the diffuse proton flux are comparable 
to observations. 



String loops that arc relatively close-by to the Milky 
Way may beam Higgses directly at the Earth and these 
would be seen as ultra-high energy cosmic rays. How- 
ever, as the particles propagate from the cusp to the 
Earth, they lose energy due to scattering with various 
components of the cosmological medium. The energy 
of a particle drops exponentially with distance from the 
cusp. If k' is the energy of a particle at the cusp, k is the 
energy at Earth, and r is the distance from the loop to 
the Earth, wc have 



k'e- r ' R 



(41) 



As the initially emitted Higgs particle decays and loses 
energy due to interactions, the products spread out over a 
wider beam angle, 9^. The total energy in all the particles 
is proportional to k9 2 and we assume that this remains 
roughly constant along the length of the beam. Therefore 
the angular spread of the beam at Earth is 



= e l h e +r ' 2R 



(42) 



where the superscripts denote Earth (E) and the loop 
(l). The number of particles at Earth with energy k 
produced from a loop of length L, at distance r follows 
from Eq. GT]) 



d$ E (fc,L,r) 



MA/ C a 4 / 3 dk' dy 1 
£2 k!V^~yTJ^yl 



dV dQ h 



(43) 



Note that the left-hand side is the flux at Earth and 
hence is at energy k, while the right-hand side contains 
the injected flux at the location of the loop and hence is 
at k'. 

The beaming solid angle is 



dn b ~ 7r(0f ) 2 = n(9 l b ) 2 e +r / R 



(/c'L) 2 / 3 



(44) 



where we have used Eq. (fl"0|) . Using Eq. ([41"]) and the 
rescalings in (|18[) and integrating over loop lengths and 
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FIG. 5: The integration region in (|45p in the ry— plane falls 
in the region bounded by the y m i n and y ma x curves from r — 
to r = T4 while the integrand depends on whether y < x or 
y > x. The dominant contribution to the integral comes from 
the region < r < n. 



spatial volume, we get 

(45) 

The limits of integration over loop lengths (y) also de- 
pend on r since (see Eq. (fT2")l ) 



2r/R 



(46) 



2/max = a mm — 
\ m' 3 



t \ = a min 



r ^/»t) (47) 



where a is defined in Eq. (|22|) . 

In the (r, y) plane of integration in Eq. (|45p. there are 
four values of r which are significant. The first, fx, is 
where y max becomes equal to a; (see Eq. (|18[1 ). and is 
given by 



e -ri/fl = 



(48) 



The second, r 2 , is where y max becomes equal to the Hub- 
ble size loop 



-r 2 /R 



k 



The third, r^, is where y m i n becomes equal to x 

e -r 3 /R = (™y /2 A 

Km) k 

and the last, r±, where y m - la equals at 



(49) 



(50) 



(51) 



The integration in the ry-plane in Eq. (|45|) can be split 
up into four integrations over the ranges (0, ri), (ri,r 2 ), 
(^2,^3) and (r3,T4). For y values less than x we can use 
yi ~ x, as argued in Sec. IIIH For y values greater than x, 
we use yt ~ y. While the full integration can be done, it 
is unnecessary because the dominant contribution comes 
from the interval (0, r\) i.e. the closest strings. The other 
integrals are suppressed by powers oik/k* (see Eq. (|29p). 
The integration over (0,ri) is 



h = 



^y ml „ y 2/i 

6 , (M\ 
In — I i? 

or \ m 



(52) 



where we have used jfj f=a x, y m ax/2/min = (M/m) 3 , and 
assumed fe/fe* <C 1. 
This leads to 



dk 



M 2 



m J (T h G^) 2 k 2 t 4 



(53) 



Using t = 10 4 Mpc = 10 17 s, and dividing by the 
Earth's surface area, 10 9 km 2 , we get the flux per unit 
area 



d<5> h 
' dAdk 



10- 



10 20 cV\ 2 /10 7 GcV\ 2 



M J 



R 



5 Mpc J 



m~ 2 s- 1 



(54) 



For M 

observed 

10- 36 m~ : 



~ 10 13 GeV, this is comparable to the 
flux of ultra high energy cosmic rays, 
s" 1 sr- 1 eV" 1 at 10 20 eV [HJ. Note that m only 
enters through the logarithm and its precise value does 
not make much difference to the overall estimate. 

The direct photon flux will also be given by Eq. (|53l) . 
but the attenuation length R will be specific to photons 
which is less than that for protons at energy 10 19 eV, and 
the branching ratio for Higgs decay into protons differs 
from that to photons. Hence the photon to proton ratio 
is 



2 
p 



N~, (R 



N v \R 



(55) 



where N 7 and N p are the number of gammas and the 
number of protons produced by a Higgs. From Figs. 9 
and 11 of Ref. and also see Fig. 2 of [H, we find 
Rj/Rp ~ 10~ 2 at 10 19 eV. If the heavy particle emitted 
by the cusp is the electroweak Higgs, the decay products 
contain a pion fraction of ~ 0.75 and a nucleon fraction 
~ 0.15. The pions then decay into gamma photons. As 
a conservative estimate of the number of photons to pro- 
tons in the decay products, we take N 1 /N p ~ 10 2 . Then 
the gamma to proton fraction at 10 19 eV is ~ 10~ 4 , 
which is consistent with the observed AUGER bound 
j/p < 0.02 [13. As pointed out in Ref. [13], at higher en- 
ergies, since R p falls quite rapidly, we expect the gamma 



8 



to proton ratio to be larger. For example, at 10 20 eV, 
Rj/Rp ~ 10 _1 and 7/p ~ 10 _1 , whereas the observed 
AGASA and Yakutsk bound is 0.3 [H]. 

We would like to point out that we have performed our 
analysis by taking R to be a constant, as in Eq. (|4Tj) . This 
is only an approximation since R depends on the energy 
of the particle and changes as the particle propagates. A 
more complete analysis would take the energy loss to be 
given by 



dk 
dr 



R(k) 



(56) 



where R(k) follows by considering the various interac- 
tions that a cosmic ray may encounter en route. Our 
results are valid only if 



dR 

~d~k 



« 



(57) 



An improved analysis should take the detailed form of 
R(k) into account. 



VII. CONCLUSIONS 

It is generally difficult to construct astrophysical sce- 
narios that can accelerate protons to high energies suf- 
ficient to arrive as ultra-high energy cosmic rays. Top 
down models involving topological defects can naturally 
give the requisite energies. Depending on the precise va- 
riety of defect, however, they are likely to suffer from an 
excess of diffuse gamma ray production. 

We have re-visited cosmic ray production from cos- 
mic strings, taking into account the possibility of a linear 
interaction between a scalar field and the string world- 
sheet. Such an interaction arises, for example, if a scalar 
field acquires a VEV ("condenses") within the string. 
Then beams of Higgs particles are emitted from cusps 
on cosmic string loops. Our analysis shows that such 
events can be responsible for the production of UHECR 
within reasonable parameters and also be consistent with 
measurements of the diffuse backgrounds and the photon 
to proton fraction in UHECR. 

Our results are also interesting because they show that 
strings with bosonic condensates are allowed in a narrow 
window of energy scales. If they are much lighter than 
~ 10 13 GeV, they will produce an excess of UHECR and 
diffuse fluxes, while if they are much heavier, they will 
cause conflicts with other cosmological constraints that 
depend on their gravitational effects. However, if the 
strings are around the Grand Unified scale, they can pro- 
duce the observed flux of ultra high energy cosmic rays 
and also be heavy enough to be detected by their gravi- 
tational signatures in the near future. 

We have focused on deriving analytical estimates of 
cosmic ray fluxes under various simplifying assumptions. 
Further work is needed to obtain more detailed estimates 
that can be compared to observations. 
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APPENDIX A: LOOP DYNAMICS AND CUSPS 

Here we summarize some known features of cosmic 
string loops that are used to obtain the radiation esti- 
mates in Sec. |TTJ 

A cosmic string loop oscillates under its own tension 
such that its world-sheet can be written as 



X(a,t) = i[ a (CT-t)+b(<T + t)] 



where a and b satisfy 



i'l = 1 = Ib'l 



(Al) 



(A2) 



where primes denote derivatives with respect to the ar- 
gument. Also, since the loop is a closed string, 



a' da = 



h'da 



(A3) 



Therefore —a' and b' arc two closed curves on a unit 
sphere whose centers of mass coincide with the center of 
the sphere. These curves generally intersect leading to 



a'(o--i) = -b'(a + t) 



(A4) 



for one or more values of a and t. Since the velocity of a 
point on the loop is 



v(tr,t) 



b'l 



(A5) 



intersection of the —a' and b' curves implies a point on 
the loop that reaches the speed of light at one instant 
per oscillation. Such a point on the string is called a 
"cusp" . The ultra-high boost factors at the cusp can be 
responsible for burst-like events from cosmic string loops. 

We choose our string paramctrization so that the cusp 
occurs at a = = t, and then expand the functions a 
and b around a cusp 



f 



1 



a(C-) = a^- + ^e + ^"C 3 



b(c + ) = b[,c + 



h'Y 2 



-W"r 3 

; D S+ 



(A6) 



(A7) 
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where £± — a ±t. The expansion coefficients are con- 
strained by 



-W ; a ' - a '>=b' -b'> = Q 



|b' | = 1 



(A8) 



These expansions give 
uj k t - k • X ~ [\/k 2 + m 2 -k)C + S kL~ 2 ( 3 + ... (A9) 



where we have taken (+ ~ £_ and denoted them collec- 
tively by £. 5 is a dimensionless constant that depends 
on the shape of the cusp. We will take S ~ 1. To have 
\oJkt— k-X| < 1, we require that both terms in Eq. (|A9[) be 
less than 1. This gives £ < min(fc/m 2 , L(fcL) -1 / 3 ). For 
k > m\/mL, we get £ < L(fcL)" 1 / 3 and for fc < mvmL, 
we get C < k/m 2 . 
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